Spin-orbit coupling and electron spin resonance theory for carbon nanotubes 
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A theoretical description of electron spin resonance (ESR) 
in ID interacting metals is given, with primary emphasis on 
carbon nanotubes. The spin-orbit coupling is derived, and the 
resulting ESR spectrum is analyzed using a low-energy field 
theory. Drastic diff'erences in the ESR spectra of single-wall 
and multi-wall nanotubes are found. For single-wall tubes, 
the predicted double peak spectrum is linked to spin-charge 
separation. For multi-wall tubes, a single narrow asymmetric 
peak is expected. 

PACS numbers: 71.10.-w, 73.63.Fg, 76.30.-v 



Electron spin resonance (ESR) serves as a valuable tool 
to experimentally probe the intrinsic spin dynamics of 
many systems. In ESR experiments one applies a static 
magnetic field and measures the absorption of radiation 
polarized perpendicular to the field direction. In the ab- 
sence of SU{2) spin symmetry breaking terms in the sys- 
tem Hamiltonian, the absorption intensity is then simply 
a (5-peak at the Zeeman energy [Q. Since spin-orbit (SO) 
interactions are generally the leading terms breaking the 
SU{2) invariance, deviations in the ESR intensity from 
the (5-peak, e.g. shifts or broadenings, are directly con- 
nected to these couplings. In this Letter we theoretically 
address the spin-orbit interaction and the resulting ESR 
spectrum for interacting ID metallic conductors, in par- 
ticular for carbon nanotubes. Nanotubes constitute a 
new class of mesoscopic quantum wires characterized by 
the interplay of strong electron-electron interactions, re- 
duced dimensionality, disorder, and unconventional spin 
dynamics |^-|6|. ESR is an important technique to elu- 
cidate aspects of this interplay inaccessible to (charge) 
transport experiments. For interacting many-body sys- 
tems, surprisingly little is known about ESR although it 
represents an interesting theoretical problem. 

Two main classes of nanotubes may be distinguished, 
namely single-wall nanotubes (SWNTs) which consist 
of just one wrapped-up graphite sheet, and multi-wall 
nanotubes (MWNTs). MWNTs contain additional inner 
shells, but transport is generally limited to the outermost 
shell Q . Evidence for the Luttinger liquid (LL) behavior 
of interacting ID electrons has been reported for charge 
transport in SWNTs ||], where one also expects to find 
spin-charge separation . Conventional wisdom holds 
that the SO coupling in ID conductors destroys spin- 
charge separation Below we show that this state- 
ment is incorrect. Indeed, the SO interaction considered 



in Ref. ^ was intended for the limited class of semicon- 
ductor quantum wires in strong Rashba and confinement 
electric fields, but in fact does not represent the generic 
SO Hamiltonian for ID conductors. The latter is derived 
below and determines the ESR intensity in SWNTs and 
MWNTs. A totally diflFerent ESR spectrum compared to 
expectations based on Ref. [Q emerges. In particular, the 
single 5-peak is split into two narrow peaks in SWNTs, 
while for the SO coupling of Ref. the spectrum forms a 
broad band with thresholds at the lower and upper edge 
[^. This qualitative difference can be traced back to 
the fact that the SO interaction in SWNTs [see Eq. (|) 
below] does not spoil spin-charge separation. Experi- 
mental observation of the peak splitting could therefore 
provide strong evidence for the elusive phenomenon of 
spin-charge separation. In MWNTs, inner shells cause a 
rather strong Rashba-type SO coupling, leading instead 
to a single narrow asymmetric ESR peak. To experimen- 
tally observe our predictions, it will be crucial to work 
with samples free of magnetic impurities whose presence 
has drastically affected previous ESR measurements for 
nanotubes [Q. 

In the standard Faraday configuration, the ESR inten- 
sity at frequency lo is proportional to the Fourier trans- 
form of the transverse spin-spin correlation function H, 



/(w) = dte''^*{S+{t)S-{0)), 



(1) 



where the static magnetic field points along the z-axis 
and S = J2i is the total spin operator. Equation (||) 
is connected to the dynamical susceptibility via I{uj) ~ 
ujx"{<^)- The Hamiltonian is H = Hq + Hz + H' , where 
Hq represents the SU (2) invariant nanotube Hamiltonian 
including electron-electron interactions, Hz = —BS^ is 
the Zeeman term and H' represents SU{2) break- 
ing terms, in particular the SO coupling. Inserting a 
complete set of eigenstates \a) oi H in Eq. (|l|), the ESR 
intensity follows as 



5(uj-{Ea-Ei,))\{a\S-\h)\ 



If H is SU{2) invariant (apart from Hz), only re- 
ceives contributions from matrix elements between eigen- 
states with equal total spin Sa = Sh. Then all states with 
= — I will contribute to form a (5-peak at frequency 
Lu = B. At zero temperature, the application of a mag- 
netic field B, taken as large enough to overcome a spin 
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gap possibly present at i? = 0, leads to a ground state 
with finite magnetization, Sq ^ 0, and the states with 
'^a — ^0 ^ ^ yield the 5-peak. Any perturbation pre- 
serving SU (2) invariance will neither shift nor broaden 
this peak, even at finite temperature, and it is therefore 
crucial to identify H'. For quantum spin chains, stag- 
gered magnetic fields and Dzyaloshinskii-Moriya interac- 
tions have been emphasized . 

Let us start with the derivation of the SO interaction 
in nanotubes. In this derivation we neglect electron- 
electron interactions which only weakly renormalize the 
SO strength 0. In a single-particle picture, the SO in- 
teraction then appears because an electron moving in the 
electrostatic potential $(r), e.g. due to the ions, sees an 
effective magnetic field v x W^. With p = mv and the 
standard Pauli matrices tr, the SO interaction reads in 
second-quantized form: 



H' = -^^ f df *t[(^x V$)-a]*. 
4m J 



(2) 



The electron spinor field ^'^(j^), defined on the wrapped 
graphite sheet, can be expressed in terms of the elec- 
tron operators Ci for honeycomb lattice site i at r^, 
= J2iXir- n) Cia, where x(^) is the 2^2 orbital 
wavefunction. The localized orbitals can be chosen as 
real-valued functions even when hybridization with 2s 
orbitals is taken into account, but their specific form is 
of no immediate interest here. We then obtain the SO 
interaction, see also Refs. [fllJll 



h.c. 



(3) 



which indeed breaks SU (2) symmetry. With obvious 
modifications, Eq. (||) applies to other ID conductors and 
thus represents a generic SO Hamiltonian. The SO vector 
Ujk = —Ukj has real-valued entries. 
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dr^{r) [Vx(r - r^) Vx(r - 4)] . (4) 



The on-site term (j ~ k) is identically zero, and since the 
overlap decreases exponentially with fj ~ fkl, we keep 
only nearest-neighbor terms in Eq. (|). We mention in 
passing that Eq. (p) has previously been found from k -p 
theory by Ando [n2|. However, his approach makes rather 
special model assumptions and is technically demanding, 
yet it does not allow to reliably compute the SO vector 
Ujk to better accuracy than specified in Eq. (Q). In addi- 
tion, the effect of SO interactions within the low-energy 
theory of nanotubes ||^ has not been analyzed. We there- 
fore take Eq. (|) as the SO Hamiltonian for SWNTs and 
MWNTs, with the SO vector (|). This formulation also 
allows to incorporate electric fields due to impurities or 
close-by electrodes in a simple and elegant manner. 

Let us first turn to SWNTs, where SO couplings are 
expected to be small. This can be rationalized from our 



approach since the SO vector (Q) vanishes by symmetry 
for an ideal 2D honeycomb lattice. A finite (nearest- 
neighbor) SO coupling can only arise due to the curva- 
ture of the wrapped sheet, stray fields from nearby gates, 
or due to defects, all of which break the high symme- 
try and in principle allow for significant SO couplings 
Focusing on the curvature-induced SO coupling for 
nonchiral SWNTs, the SO vector only depends on bond 
direction, = Ua, for the nearest-neighbor bonds 

6a (a ~ 1, 2, 3) of the graphite sheet ^. 

To make progress, we employ the effective field the- 
ory approach ||]. Neglecting the (here inessential) "fla- 
vor" index due to the presence of two Fermi (K) points 
Ho then corresponds to a spin- 1/2 LL described by 
charge/spin interaction parameters Kc < l,Ks and ve- 
locities Vf.^^ — vp/K^/s with the Fermi velocity vp = 
8 X 10^ m/sec; SU{2) invariance of Hq fixes Kg — 1 
[ |l3t . The LL Hamiltonian completely decouples when ex- 
pressed in terms of spin and charge boson fields [ p^ , and 
the Zeeman term Hz only affects the spin sector. Writ- 
ten in terms of right- and left-moving fermions ippi/L, the 
continuum version of Eq. (^) is (up to irrelevant terms) 
H' = Hi+ H2 with 



Hi = j dxX-{JL' Jb), 

H2 ~ dx ^ tplX'-aidxipr + ^-C- 

•' r=R/L 

With the unit vector et along the tube axis, we use 



(5) 
(6) 



A = 2 Im^e' 



(et • Sa)Ua 



and neglect oscillatory terms which average out on large 
length scales. These oscillations are governed by the 
wavevector 2kp corresponding to the doping level /i, 
kp = \^\/vp, with typical values |^| « 0.3 to 0.5 eV 
0. Finally, Jr^l — ''P'^r/l^^ /'^)'4^p!^/l standard 
SU{2) spin currents. The perturbation (||) has scaling 
dimension 1 (relevant) whereas Eq. (^) has dimension 2 
(marginal). Therefore the leading SO contribution re- 
tained in what follows is Eq. (||). 

Remarkably, the SO interaction (||) acts exclusively in 
the spin sector and hence does not spoil spin-charge sepa- 
ration. As a consequence, since electron-electron interac- 
tions affect only the charge sector, the ESR intensity can 
be computed using the equivalent fermionic spin Hamil- 
tonian 



Hf = 



dx 



-irvsipld^ipr + K ■ Jr 



(7) 



where A± = -B ± A. Since Hf is bilinear in the fermions, 
after some straightforward algebra, the exact ESR spec- 
trum follows for arbitrary temperature 
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with A± = |A±|. As a consequence of SO couphng, the 
single i5-peak splits into two peaks but there is no broad- 
ening. The peak separation is |A+ — A_|, and the peak 
heights are generaUy different. To lowest order in X/B, 
the two peaks are located symmetrically around uj — B. 
Notice that for B perpendicular to the effective SO vector 
A, the splitting is zero. It should be stressed that these 
results hold both for the non-interacting and the inter- 
acting case. The double peak is therefore not directly re- 
lated to the doubling of singularities in the single-electron 
Greens function commonly associated with spin-charge 
separation. However, for the interacting case realized in 
SWNTs ||] , the double peak structure is only possible if 
spin-charge separation is present. Otherwise the charge 
sector will mix in, leading to broad bands with threshold 
behaviors Closer inspection shows that inclusion of 
the subleading term preserves the splitting into two 
peaks, but the peaks now acquire a finite width ~ |A'|. 

Experimental observation of the predicted double-peak 
spectrum could then provide strong evidence for spin- 
charge separation. In practice, to get measurable inten- 
sities, one will have to work with an ensemble of SWNTs. 
The proposed experiment may be possible using electric- 
field-aligned SWNTs, or by employing arrays of identical 
SWNTS ||l5|. In more conventional samples containing 
many SWNTs, however, the SO vector A can take any 
direction. Assuming a uniform probability distribution 
for the orientation of A, the average can easily be done. 
For T = and A' = 0, we find 



[{uj + Bf - X 
IGvsB^X 
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+ X-LU)0{UJ- \B-X\) 



(9) 



with the Heaviside function 9{z). For B > X, this asym- 
metric spectrum has the width Aw = 2A, which in turn 
allows to extract the SO coupling strength A from ESR 
measurements. 

For the remainder, we then focus on MWNTs where 
we first contemplate a simple two-shell model. Exper- 
imental evidence ^] is consistent with the assumption 
that inter-shell tunneling is strongly suppressed. There- 
fore doping due to charge transfer from the substrate, 
the attached leads, or due to oxygen absorption, should 
only affect the outermost shell, while « for inner 
shells. Under this assumption, basically all conduction 
electrons contributing to the ESR signal reside in the 
outermost shell. Moreover, the electrostatic potential $o 
of the outer shell differs from the inner-shell potential 
$j. In effect, we can then restrict attention to the out- 
ermost shell (of radius R) alone, but in a radial electric 
field of size E « 2ci2A$/i?, where Ci2 is the inter-shell 



capacitance per length and A$ = $i — $o. The general 
expression (||) for the SO vector then yields after some 
algebra for a given bond Sa 



Ua = {u/vf)E X (5q, 



(10) 



where E is a. unit vector perpendicular to the tube sur- 
face, and u K, ci2eA$/(m^(ii?) with the C-C distance 
d= 1.42 A. 

To proceed, we turn to the low-energy theory, again 
for only one K point. The influence of interactions 
is expected to be less dramatic in MWNTs compared 
to SWNTs, and here we focus on the non-interacting 
case. In real space, we then have a Dirac Hamiltonian, 
Hq = —ivp Jdr'i>''{T- V)^', where the integral extends 
over the tube surface and the Pauli matrices r act in 
sublattice space. The SO contribution appropriate for 
MWNTs follows by inserting Eq. (|l^) into Eq. ^), which 
yields the manifestly Hermitian term 



H' 



-lUVp 



df vI'Vr o-^ 



(11) 



We take the magnetic field parallel to the tube axis and 
include both orbital and Zeeman contributions . 

The full Hamiltonian can then be diagonalized. The 
dispersion relation contains four branches, e(/c) — 
±e±(k), where the ± signs are independent and 



e±(fc) 



vj.{k^ + Q^) + {B/2)^±VF 



(12) 



1/2- 
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with Q = 2^/2u measuring the SO strength and k = \k\. 
Using Eq. (|l2|), the ESR intensity at finite temperature 
T reads 



-1 



1 - e-"/^ / Stt 



dk H(e+, A:; uj) S{uj — e+ + e_) 



-(4 



c2 ^2 



[nF(e+) - n_F(-e+) - npie^) + nF(-e-)] 



(13) 



where np{e) = 1/(1 + e^'^^*^^/-^) is the Fermi-Dirac dis- 
tribution function. The integral in Eq. (|l^) includes an 
integration over the momentum parallel to the MWNT 
axis and a discrete summation over the quantized trans- 
verse momenta, k± = {n — 4>)/R, where (p — nR^B/{h/e) 
is due to the orbital effect of the applied magnetic field 
and the summation extends over integer n. Furthermore, 



with u' = Lo' -u and /(w', k) = v^-k"^ - (-S/2 ^ 

H(c^',fc;co) = 2/(L^',fc)/(-cS',fc)x 
[{B/2 - J){B/2 + u') - vlk^] 
~2vlQ^ [{B/2 - u'){B/2 - u)')f{uj\ k) 
+ {B/2 + uj'){B/2 + Lo')f{-uj', k)] . 
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As a simple check, for Q 



(14) 

one recovers the expected 



(5-peak from Eq. (13) 
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The result (|T|) can be understood in simple physical 
terms. The ESR intensity receives contributions from 
transitions between states of energy e_ to states of energy 
e+, and each contribution is weighted by a factor which 
takes into account the occupation of the levels. In order 
to arrive at Eq. (|l3|), we have neglected a contribution 
coming from transitions between — e± and e=p states. This 
is consistent because these transitions only contribute for 
very large frequencies to > while the frequency scales 
relevant to ESR are much lower. In addition, the signal 
given by this contribution is very small in comparison to 
the term that we keep. 

Inspection of Eq. shows that the ESR spectrum 
of a MWNT at low temperatures contains only a sin- 
gle narrow asymmetric peak, whereas more structure ap- 
pears at higher temperature due to the activation of the 
transverse subbands [0. Here we focus on the low- 
temperature ESR spectrum, which is shown in Fig. |l| for 
typical parameters. The peak has an asymmetric line- 
shape which strongly depends on temperature. At zero 
temperature, the intensity maximum is at the frequency 



Wo = ez 



1 



-0{{ez/t^f 



(15) 



where ez — \/ + 2vpQ^. As is apparent in Fig. |^, 
when increasing the temperature, the position of the 
maximum slowly moves to smaller frequencies. The 



linewidth is of the order of 



(16) 



To conclude, we have presented a theoretical descrip- 
tion of the spin-orbit coupling and the ESR spectrum for 
ID conductors, in particular for carbon nanotubes. In 
SWNTs, spin-charge separation should not be affected 
by spin-orbit coupling, and hence we expect a double 
peak. The peak distance and their height provides infor- 
mation about the spin-orbit coupling strength, and their 
width points to violations of spin-charge separation. The 
ESR spectrum of a MWNT only exhibits a single narrow 
peak, whose location, lineshape and linewidth provide 
information about the Rashba-type spin-orbit coupling 
and intrinsic electric fields. The ESR spectra of SWNTs 
and MWNTs are therefore fundamentally different and 
reflect distinct mechanisms of spin-orbit coupling. 
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FIG. 1. Typical ESR intensity of a MWNT at low temper- 
atures. Parameters in the plot are </> = 0, p = 0.1, B — 0.0014 
and vfQ = 0.05, where energies are measured in units of 
2nhvF /R. The value of B corresponds to a field of 10 Tesla. 
Note the frequency units, pointing to a very narrow ESR peak. 
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